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Abstract 

We consider the behaviour of the perturbative QCD corrections to the 
Re+e- ratio, in the hmit that the cm. energy ^/s vanishes. Writing Re+e- (s) = 
3^^Qj(l + lZ{s)), with Qf denoting the electric charge of quark flavour /, 
we find that for Nf < 9 flavours of massless quarks, the perturbative cor- 
rection TZ{s) to the parton model result smoothly approaches from below 
the infra-red limit 7?.(0) = 2/6, as s— i>0. Here h = (33 — 2Nf)/6 is the first 
QCD beta-function coefficient. This freezing holds to all-orders in perturba- 
tion theory. The s-dependence can be written analytically in closed form in 
terms of the Lambert W function. 
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In the ultra-violet s^oo limit of QCD the renormalized coupling vanishes, 
and this property of Asymptotic Freedom underwrites the successful use of 
perturbative methods in testing the theory In the infra-red limit s—>-0, 
however, one may expect that perturbation theory will break down , with 
typically a Landau pole singularity in the coupling when s~Aq^^, and that 
non-perturbative effects will be important. However, the phenomenological 
virtues of assuming a frozen couplant, with the renormalized as{s) approach- 
ing a constant value as(0)/7r ~ 0.3 in the infra-red have long been recognised 
[2-5]. In a pioneering paper Mattingly and Stevenson investigated the be- 
haviour of the perturbative corrections to Re+e- including third-order QCD 
corrections, in the framework of the Principle of Minimal Sensitivity (PMS) 
approach [0]. Their PMS optimized coupling indeed froze to a value around 
0.26 below 300 MeV. These predictions were then smeared using the tech- 
nique of Poggio-Quinn- Weinberg (PQW) 0, and were in suprisingly good 
agreement with similarly smeared experimental data for Re+e-- Some scep- 
ticism about the existence of infra-red fixed point behaviour had previously 
been expressed ||^. In this letter we wish to demonstrate that including all- 
orders in perturbation theory the perturbative corrections to Re+e- do freeze 
in the infra-red. The limiting value being 2/6, where b = (33 — 2Nf)/6 is the 
first beta-function coefficient of QCD with Nf quark flavours. We assumed 
massless quarks and for freezing to this limit one requires Nf < 9 flavours. In 
fact the freezing behaviour does not correspond to an infra-red flxed point in 
the beta-function, but rather stems from the energy dependence induced by 
analytical continuation from the Euclidean to Minkowskian region in defln- 

ing Re+e-- 

We begin by deflning the Re+e- ratio at cm. energy a/s, 

Re+As)=^^^^^^^j^^ = 3^g/i + n{s)) . (i) 

Here the Qf denote the electric charges of the different flavours of quarks, 
and 7^(s) denotes the perturbative corrections to the parton model result, 
and has a perturbation series of the form, 

7^(s) = a + ^r,a"+^ . (2) 

Here a=as{fi^) /'n' is the renormalized coupling, and the coefficients vi and r2 
have been computed in the MS scheme with renormalization scale fi"^ = s 



2 



0, Re+e~ is directly related to the transverse part of the correlator of two 
vector currents in the Euclidean region, 

{Qt^qu - g^uQ^ms) = ^Ti^t j ciVe'"-" < 0|T[j^(x)j,(0)]|0 > , (3) 

where s = — > 0. To avoid an unspecified constant it is convenient to take 
a logarithmic derivative with respect to s and define the Adler D-function, 

D{s) = -s-^m ■ (4) 

This can be represented by Eq.(l) with the perturbative corrections TZ{s) 
replaced by 

V{s) = a + ^c/na"+^ . (5) 

n>0 

The Minkowskian observable TZ{s) is related to I){—s) by analytical contin- 
uation from Euclidean to Minkowskian. This can be elegantly formulated 
as an integration around a circular contour in the complex energy-squared 
s-plane ilig, 

7^(s) = — / devise''') . (6) 

Expanding J){se^^) as a power series in a=a(se*^), and performing the 6 in- 
tegration term-by-term, leads to a "contour-improved" perturbation series, 
in which at each order an infinite subset of analytical continuation terms 
present in the conventional perturbation series of Eq.(2) are resummed. It 
is this complete analytical continuation that builds our claimed freezing of 
Tl{s). In contrast Ref.0] used the conventional fixed-order perturbative ex- 
pansion of Eq. (2) in the PMS approach. It is useful to begin by considering 
the "contour-improved" series for the simplified case of a one-loop coupling. 
The one-loop coupling will be given by 

a(s) = ^ . (7) 

As described above one can then obtain the "contour-improved" perturbation 
series for 7l{s), 

oo 

n{s) = Ai{s) + J2dnAn+i{s) , (8) 

n=l 
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where the functions are defined by, 



2ttJ_^ 27rJ_^ [1 + z66'a(s)/2J 



(9) 



This is an elementary integral which can be evaluated in closed- form as 11 
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We then obtain the one-loop "contour-improved" series for TZ{s] 
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[l + bH^a^{s)/Ay 
(11) 

The first arctan term is well-known, and corresponds to resumming the in- 
finite subset of analytical continuation terms in the standard perturbation 
series of Eq.(2) which are independent of the dn coefficients. Subsequent 
terms corrrespond to resumming to all-orders the infinite subset of terms 
in Eq.(2) proportional to di,d2, . ■ ., etc. It is crucial to notice that in each 
case the resummation is convergent, provided that \a{s)\ < 2 /'Kb. Thus, even 
though the series in Eqs.(2,5) are divergent because of the n! growth of the 
dn coefficients, the resummations implicit in performing the integration of 
Eq.(6) term- by-term are well-defined. In the ultra-violet s oo limit the 
An{s) vanish as required by asymptotic freedom. In the infra-red s — >■ limit, 
the one-loop coupling a{s) has a "Landau" singularity at s = A|j^. However, 
the functions A„(s) resulting from resummation, if analytically continued, are 
well-defined for all real values of s. Ai{s) smoothly approaches from below 
the asymptotic infra-red value 2/6, whilst for n > 1 the An{s) vanish. Thus, 
as claimed, 7^(s) is asymptotic to 2/6 to all-orders in perturbation theory. It 
might be thought that the existence of such an infra-red limit, independent 
of the higher-order structure of perturbation theory, is remarkable. In fact, 
however, it is to be expected. A useful analogy is the exponentiation to all- 
orders of large infra-red logarithms which appear for jet observables such as 
thrust distributions in e"'"e~ annihilation Here the standard fixed-order 
perturbation theory breaks down in the two-jet region as these logarithms 



+. 
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become infinite. However, if tlie large logaritfims are resummed to all-orders 
one builds an exponential factor , and the thrust distribution smoothly ap- 
proaches zero in the two-jet region, to all-orders in perturbation theory. We 
should note a subtlety in our derivation above. As defined by the integration 
in 9 around the circle in Eq.(9) the An{s) are not defined for s < A|j^ . 
However the v4„(s) given by Eq.(lO) are defined for all real s. A more careful 
derivation would use instead of the contour integral as a starting point, the 
dispersion relation, 

which clearly avoids the "Landau" singularity and is well-defined for all real 
s. One can directly obtain the result of Eq.(lO) from the dispersion relation 
by simple manipulations. It should be noted that this route is equivalent to 
the Analytic Perturbation Theory (APT) approach, in which it is advocated 
that Minkowskian observables are expanded in a basis of functions 2l„(s) ob- 
tained by performing an integral transform of the Euclidean coupling using 
the dispersion relation of Eq.(12). In the case of the Re+e- ratio the resulting 
APT expansion corresponds to the "contour-improved" expansion of Eq. (8) 
and 2l„(s) = A„(s), given by Eq.(lO). The infra-red freezing and absence 
of a "Landau pole" in the 2l„(s) has previously been discussed in the APT 
approach, and provides one of its major motivations. Analytic expressions 
for the one-loop 2l„(s) have been given, and freezing to all-loops in APT can 



be demonstrated. For a recent review see Ref. 121, and references therein 



We now move to the more challenging problem of what happens for real- 
istic QCD beyond the simple one-loop approximation. The freezing is most 
easily analysed using a renormalization scheme in which the beta-function 
equation has its two-loop form, 

|^ = 4^V)(l + ca(/.^)). (13) 



This corresponds to a so-called 't Hooft scheme in which the non- 
universal beta-function coefficients are all zero. Here c = (153 — 19Nf)/12b 
is the second universal beta-function coefficient. The key feature of these 
schemes is that the coupling can be expressed analytically in closed-form in 
terms of the Lambert W function , defined implicitly by W {z)ex]i{W (z)) = z 
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TM. One has 



/i 



-b/c 



(14) 



where Ajj^ is defined according to the convention of ||T5[ , and is related to 
the standard definition |jl6[ by Aj^g = (2c/6) '^^''Ajg^. The "—1" subscript on 
W denotes the branch of the Lambert W function required for Asymptotic 
Freedom, the nomenclature being that of Ref. [1^. Assuming a choice of 
renormalization scale /i^ = , where a: is a dimensionless constant, for the 
perturbation series of T'(s) in Eq.(5), one can then expand the integrand in 
Eq.(6) for TZ{s) in powers of a = a(xse*^) , which can be expressed in terms 
of the Lambert W function using Eq.(14), 



-1 



c[l + W{A{s)eiK^)] 



where 
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The functions v4„(s) in the "contour-improved" series are then given, using 
Eqs(15,16), by 



An{s) = — 
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Here the appropriate branches of the W function are used in the two regions 
of integration. As discussed in Refs.[|T^, by making the change of variable 
w = W{A{s)e^^^) we can then obtain 

(18) 
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This is an elementary integral and noting that Wi{A{s)e~^^'^) = [W^i{A{s)e'^^'^ 
we obtain for n = 1, 
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(19) 
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where the 2/6 term is the residue of the pole at w = 0. In Ref. [T^l this 
contribution was omitted in error. For n > 1 we obtain 



c^/Ttt [ \l + W-i{A{s)e^-) J frik{l + W_i{A{s)e'^^))'' 

(20) 

Here the contributions from the poles at w = and w = —1 cancel exactly. 
Equivalent expressions for the have been obtained in the APT approach 

[ll9| . Provided that b/c > 0, which will be true for Nf < 9, the functions 
An{s) are well-defined for all real values of s. As s ^ oo they vanish as 
required by Asymptotic Freedom. As s — > 0, Ai{s) smoothly approaches 
the infra-red limit 2/6 from below, as the second term in Eq.(19) vanishes 
in the limit. Whilst for n > 1 the A„(s) vanish as s ^ 0, and so, as in the 
one-loop case, 7?.(s) is asymptotic to 2/6 to all-orders in perturbation theory. 
The cancellation of pole contributions noted above is crucial in achieving 
this. We should point out that for 6/c > 4 we move to other branches of the 
Lambert W function in order to keep a(/i^) continuous. This just changes 
the value of the branch of the Lambert W function in A„(s), and will not 
alter our result for s — >■ 0. We find b/c> 4 first occurring for Nf = 7. Again 
we need to refine the above argument slightly, using the dispersion relation 
of Eq.(12), which after a change of variable w = W{A{t)) yields the integral 
of Eq.(18). As in the one-loop case the of Eqs. (19,20) can be obtained 

as all-orders resummations of subsets of terms in the standard perturbation 
series of Eq.(2), these resummations are convergent provided \a{s)\ < 2/7r6. 
In Figures 1-3 we plot the functions Ai{s),A2{s) and ^3(5), respectively, as 
functions of sx/ A^j^. Nf = 2 flavours of quark are assumed. 



7 



0.0004 



0.0002 



-0.0002 



-0.001 




Figure 3: As Fig.l but for ^(s) of Eq.(20). 



We have shown that the freezing occurs to all-orders in perturbation 
theory, and thus must occur independent of the choice of renormalization 
scheme. The use of a 't Hooft scheme turns out to make the freezing man- 
ifest. We used a general choice of renormalization scale /i^ = xs. The true 
infra-red s-dependence of 7^(s) does not depend on the unphysical parame- 
ter X. One should rather eliminate the /i-dependence of the result altogether 
by completely resumming all the ultra-violet logarithms which build the s- 
dependence. This so-called Complete Renormalization Group Improvement 
(CORGI) approach |^ corresponds to choosing fx^ = e~'^'^/^s, where d is the 
NLO perturbative correction di to T>{s) in Eq.(5), in the MS scheme with 
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jj^ = s. One then has the "contour- improved" CORGI series, 



n{s) = A^{S) + Y,XnAn+l{s) , (21) 



n=2 



where the X„ are the CORGI invariants, and only X2 is known. Now A{s) = 
(— l/e)(A/i/A£)) where AD=e'^^^Ajjg. For an infra-red fixed point at 
7?.(0) = 71*, corresponding to a zero of the beta-function, one expects the 
asymptotic behaviour [Q] 

7^(s) - 7^*~s^ , (22) 

where 7 is a critical exponent. Our freezing to 2/6 is instead driven by the 
analytical continuation from the Euclidean to Minkowskian regions, and one 
obtains the asymptotic behaviour, 

, , 2 -l/c-2/6 , , 



Again this involves the ubiquitous Lambert W function. 

We should stress that, of course, the result TZ{0) = 2/b is in itself only 
of abstract interest since s = 4m^ is the threshold in full QCD with mas- 
sive quarks. The important conclusion is that while the conventional per- 
turbation series of Eq.(2) breaks down at the spurious Landau pole in the 
coupling a{s), this is eliminated by completely resumming all the analytical 
continuation terms, so that the "contour-improved" (or APT) perturbation 
series is well-behaved in the infra-red. It was crucial in investigating this to 
be able to evaluate the functions v4„(s) in closed analytic form. In previous 
phenomenological investigations [T^ these functions were evaluated by nu- 
merical integration with Simpson's Rule, making it impossible to go further 
into the infra-red than the Landau pole obstruction. Whilst the infra-red 
freezing is most easily investigated using the "contour improved" or Ana- 
lytic version of the perturbation series in Eqs.(8,21), it should be stressed 
that the freezing is an all-orders result of QCD perturbation theory. The 
standard expectation is that by itself the all-orders perturbation series is 
ambiguous due to the presence of infra-red renormalons, these renormalon 
ambiguities cancelling against corresponding non-logarithmic UV divergences 
in the non-perturbative Operator Product Expansion (OPE) [^. The van- 
ishing of the An{s) for n > 1 in the infra-red means that these renormalon 
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ambiguities also vanish, and so the imphcation is that the non-logarithmic 
UV divergences of the OPE also vanish in the infra-red. This of course says 
nothing about the infra-red limit of the resummed OPE, indeed we know that 
there are important non-perturbative effects which build a complicated set 
of hadronic resonances, but the interesting observation is that perturbative 
and non-perturbative effects are separately well-defined in the infra-red limit 
for Minkowskian quantities. This is not the case for Euclidean quantities 
where the Euclidean APT coupling necessarily includes a resummation of 
non-perturbative OPE terms . 

There are evidently many phenomenological applications of the "contour- 
improved" CORGI perturbation series for TZ{s), in particular one can repeat 
the analysis of Ref. and compare PQW smeared data for i?e+e-("S) with 
the similarly smeared perturbative freezing. This exercise has been performed 



in the APT approach in Ref.|2^, and good agreement found. There will also 
be applications to estimating uncertainties in a{Mz), and in estimating the 
hadronic corrections to the anomalous magnetic moment of the muon. We 
hope to report on these aspects in a future publication. 
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